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We study the advantage of the co-existence of future ground and space based gravitational wave 
detectors, in estimating the parameters of a binary coalescence. Using the post-Newtonian waveform 
for the inspiral of non-spinning neutron star-black hole pairs in circular orbits, we study how the 
estimates for chirp mass, symmetric mass ratio, and time and phase at coalescence are improved by 
combining the data from different space-ground detector pairs. Since the gravitational waves pro¬ 
duced by binary coalescence also provide a suitable domain where we can study strong field gravity, 
we also study the deviations from general relativity using the parameterized post-Einsteinian frame¬ 
work. As an example, focusing on the Einstein telescope and DECIGO pair, we demonstrate that 
there exists a sweet spot range of sensitivity in the pre-DECIGO phase where the best enhance¬ 
ment due to the synergy effect can be obtained for the estimates of the post-Newtonian waveform 
parameters as well as the modification parameters to general relativity. 
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I. INTRODUCTION 

General relativity (GR) has enjoyed a very successful run as a theory describing the classical Universe by passing 
all the tests available so far [1]. Yet two of the unique predictions of GR, the black holes (BHs) and gravitational 
waves (GWs), are waiting for direct confirmation. There is an ever increasing indirect evidence of stellar mass BHs as 
well as super massive BHs at the center of each galaxy [2]. Similarly, the Universe is immersed in a GW background, 
analogous to the cosmic microwave background, but we still have no glimpse of it. Thanks to several decades of 
research developing detectors with ever increasing sensitivity, we are likely to be very close to detecting GWs and to 
open a new era in astronomy. It will not be an overstatement to say that if the 20th century belonged to astronomy 
with electromagnetic waves, the 21st century would definitely belong to astronomy with GWs [3, 4]. The detection 
of GWs would not only be a triumph of GR as the theory describing our classical world, but it would also give us a 
unique opportunity to detect BHs by directly mapping their spacetime. 

GR has been subjected to a plethora of experimental tests, but mostly the explored regimes are those in which the 
gravitational field is weak and particle velocities are small relative to the speed of light [1]. Thus, one of the most 
promising outcomes of GW astronomy would be the possibility of testing the validity of GR in the dynamical strong 
field regime [5] . The prospect of discovering that GR may need some modification in the strong field regime has also 
motivated the development of several alternative theories of gravity. GW signals will provide a testing ground for 
these theories and a wide range of tests of GR have already been proposed. Tests that use GWs coming from compact 
binary coalescence include those proposed in [6-12]. The coalescing binaries composed of neutron stars (NS) and BHs 
are amongst the most likely sources for GW signals to be observed by ground-based interferometric detectors. When 
observed, these GW signals will probe the strong field regime of gravity and allow the implementation of these tests 
of GR. 

The post-Newtonian formalism (PN) has been used to model the inspiral part of the binary evolution within GR. 
In the PN formalism physical quantities of interest such as the conserved energy, flux etc. are found as expansions in 
the small parameter (v/c), where v is the characteristic speed of the binary system and c is the speed of light [13]. In 
standard convention 0((v/c)") corrections counting from the leading order are referred to as a (n/2)PN order term. 
The expansion of the GW waveform is currently known up to 3.5PN, and the binary motion is known up to 4PN order 
[13]. Inspiraling binaries are suitable for studying the strong field regime. The orbital velocities in these systems can 
go as high as 0.5c, and therefore higher order PN corrections are relevant. The standard data analysis techniques used 
for the detection and analysis of GW signals depend on the availability of accurate templates to identify weak signals 
buried in the noise. However, in order to carry out the tests of GR, it is impractical to make accurate templates for all 
possible alternative theories of GR. A more feasible way of carrying such tests is to adopt general non-GR templates 
to model the signals. Arun et ah, proposed such templates in which the expansion coefficients of the phase in the 
Fourier transform of the GW waveform (see equation (3) ahead) are treated as fitting parameters [9]. In GR each 
coefficient is specified by the masses and spins of the constituting compact objects of the binary. In principle, this 
relation can be completely different in some alternative theory and may even involve other parameters. 
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Yunes and Pretorius developed the parameterized post-Einsteinian (ppE) framework, which allows for character¬ 
izing a wider range of deviations to the amplitude and phase of the waveform [10]. Similar to the parameterized 
post-Newtonian framework they introduced ppE parameters, but instead of parametrizing the metric tensor, they 
parametrize the GW response function using a generic template family (see section IIB ahead). This family can 
accommodate the inspiral phase of most of the known alternative theories with appropriate choices of the parameters. 
Cornish et al. applied the ppE approach to simulated data to determine the level at which departures from GR can 
be detected, and also analysed the bias introduced in the extraction of the model parameters due to the assumption 
of the wrong theory [11]. 

In this work we study the synergy between the ground and space based GW detectors in estimating the parameters 
of the inspiraling binaries, to see whether there is any gain in combining measurements from two detectors. For this 
we use the ppE framework developed by Yunes and Pretorius [10]. We consider two future space based detectors 
(DECIGO and eLISA) and two future ground based detectors (Advanced LIGO and ET). As representative test cases 
we consider two binary systems, a 1.4 Mq -|- 10 Mq binary and a 1.4 Mq + 100 Mq binary, where Mq is the solar 
mass. Within the PN framework, we study the synergy effect in detail for the chirp mass, the symmetric mass ratio 
and the time and phase at coalescence, and within the ppE framework we study the synergy effect on the estimation 
of the ppE parameter characterizing deviation from GR. The plan of the paper is as follows: in the next section we 
briefly describe the waveforms used in this analysis (II A, IIB), the noise curves for the different detectors (IIG), and 
the methodology for error estimation (IID). We summarise our results in section III. 


II. METHODOLOGY 


A. PN formalism 


We use the PN expression for the Fourier transform H{f), of the GW signal coming from the inspiral of binary 
compact objects in a circular orbit (under the stationary phase approximation) [14]. Moreover, the restricted PN 
waveforms are used which keep the higher terms in phase but only take the leading terms for the amplitude [15]. 
This is because in the matched filtering analysis used for the parameter estimation of binaries the correlation of two 
waveforms is more sensitive to deviation in phase than the amplitude. This is not true for binaries with misaligned 
spins, where the amplitude modulation on precession time scale is also important to determine the spin parameters 
[16]. As we consider non-spinning binaries for simplicity, this modulation is beyond the scope of this paper. 

Hence we consider only the leading order term in the amplitude, while the phase terms are taken up to 3.5PN order 
(note that terms with exponent k = n are {nl2) PN order terms). Namely, we adopt the waveform 

H{f) = Af~'^^^ exp(iV'(/) + *7r/4), (1) 


where the Fourier amplitude A is given by 


and the phase if{f) is given by 


A = 


C 

Dtt^A 



( 2 ) 


'*/'(/) = 27r/tc + (l)c + 
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( 3 ) 


Here C is an 0(1) dimensionless geometric factor that depends on the relative orientation of the binary and the 
detector (average over all orientations C = 2/5), v is the symmetric mass ratio, M is the total mass of the binary, 
D is the luminosity distance and tc and (f>c are the time and phase at coalescence, respectively. For the purpose of 
this paper, the explicit expression for C is not necessary as we also neglect the effects of the orbital motion of the 
space antenna. The expressions for the amplitude and the phase are often also expressed in terms of the chirp mass 
M = as 


A = 
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The values of the coefficients in the PN expansion of the Fourier phase are as follows [14]: 
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As mentioned earlier, we consider two binary systems, a 1.4 Mq + 10 Mq NS-BH binary and a 1.4 Mq + 100 Mq 
NS-BH binary. 


B. Corrections due to modified gravity: ppE formalism 


One of the most important outcomes of GW detection would be placing bounds on alternative theories of gravity [17]. 
Several such theories have been proposed to explain the observed late time acceleration of the Universe. Within the 
framework of GR the explanation is given by assuming exotic matter fields, but if classical GR is realised as a low 
energy limit of some more fundamental theory, we might expect non-trivial modifications to it {e.g. non-minimally 
coupled dilatons etc.). 

GR has been tested with high accuracy in solar system and also through binary pulsars, and in both cases grav¬ 
itational fields are weak and particles velocities are small [1]. The lack of experimental verification of GR in strong 
fields may not be a problem for long once we start to observe GWs directly. In particular, GW detectors in space like 
eLISA and DECIGO will be able to map the spacetime geometry of super-massive BHs at the center of galaxies. This 
can be done very accurately by observing GWs emitted by their inspiraling satellites, allowing us to test alternative 
theories of gravity. 

For studying the modification to GR using a parameterized approach we follow the ppE formalism. The modification 
to the Fourier transform H{f) of the GW signal (under the stationary phase approximation) from the inspiraling 
binary (1) would be characterized by the leading order corrections to the Fourier amplitude and the phase as 


A = (^1 + ^ AgrU), (7) 

and 

V’(/) = ■ (8) 

Here AoRif) is given by equation (2), and ipGRif) is given by equation (3), the standard GR expressions. The 
coefficients g^ and /3i, in general, may depend on parameters such as chirp mass, spin angular momentum etc. As 
explained in previous subsection, for the present analysis we concentrate on the modification in the phase only, 
meaning A{f) = AcRif), i.e., g* = 0. 


C. Noise curves 

The output of a GW detector s{t), can be decomposed into two components: the GW signal h{t), and the detector 
noise n(t), i.e. s{t) = h{t) + n(t). For simplicity, we assume that the detector noise is stationary and Gaussian. 
Stationarity implies that the different Fourier components of the noise are uncorrelated, and hence the one-sided noise 
power spectral density (PSD) «S'„(/) can be defined by 
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where n{f) is the Fourier transform of the noise n(t). A commonly used quantity to describe sensitivity curve is 
square root of PSD in units of and when integrated over positive frequencies it gives mean square amplitude 

of the signal in detector [18]. We describe the noise spectral densities of all the detectors used in this study below. 


Future space based detectors 


DECIGO: The Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is a future plan of a space 
mission for observing GWs around / ^ 0.1 — lOHz, proposed by [19]. The following noise curve is taken from [20] 


Sn{f) = 6.53 X 10-4® 



+ 4.45 X 10-®4 


('AV* 1 


+ 4.94 X 10-^2 



(9) 


where fp = 7.36Hz. 

eLISA; eLISA is a slightly down-graded mission which is derived from the Laser Interferometer Space Antenna 
(LISA) proposal. This mission is targeted at the low-frequency band: O.lmHz to IHz. The following noise curve is 
taken from [21]: 


with 


^ 20 4^„ec(/)/(2^/)4 + S^nif) + Somnif) f f V 

4 L2 0.41c/2L J 


Saccif) = 2.13 X 10-29(1 + 10-4Hz//)m2s-4Hz-4, 
= 5.25 X 10-2®m2Hz-4, 

Somnif) = 6.28 X 10-23m2Hz-4. 

where L = 10®km and c is the speed of light. 


( 10 ) 


( 11 ) 


Future ground based detectors 

Advanced LIGO: As the representative of the second generation of ground-based Laser Interferometer 
Gravitational-Wave Observatory (LIGO) The noise curve for the advanced LIGO detectors is taken from [22]. 


SM) = 10 -"® 



( 12 ) 


where /o = 215Hz. 

ET: The Einstein Telescope (ET) is a European- commission project which aims at developing the third generation 
Gravitational Wave observatory. The following noise curve is taken from [22]: 
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where /o = lOOHz. The noise curves for the various detectors are plotted in figure (1). 
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FIG. 1: Noise spectra for the different GW detectors used in this work. 


D. Likelihood function and Fisher matrix 


In this section we briefly overview parameter estimation and Fisher matrix method of error estimation. For a review 
on this topic, we refer to [23]. We assume that the GW signal depends on the parameter vector 6. In case of the PN 
waveform 9 — {M, </>c}, and in case of the ppE framework 9 — {A4, v, tc, (/'c, A} (in this work the ppE parameter 

at is fixed to zero for all i). Now the noise weighted inner product for two signals hi(t) and /i 2 (t) is defined as: 


(^1A 2 ) — 2 



kif)h2if) +Kif)hi{f) 

Sn{f) 


df, 


(14) 


where hi{f) and Ii 2 (/) are the Fourier transforms of hi{t) and /i 2 (t), respectively, and * represents the complex 
conjugation. The inner product is defined so that the probability that the signal depends on the parameters 9 
(assuming the noise is Gaussian) is given by 


P{s\9) oc 

Hence, the inner product is similar to defining a chi-squared merit function. Given an output s(t), the best fit GW 
waveform h{9) is obtained by minimizing this inner product ([23]). For different realizations of noise, we may obtain 
slightly different values of the parameters, but for large signal to noise ratio they will all be centred around the correct 
values say 9 with some spread A0. The estimation errors follow a Gaussian distribution 

P(A6»0 cx (^ 5 ) 

where F^ is the Fisher matrix defined as ([15]) 


“ Ua dej- 


(16) 


The covariance matrix C is the inverse of the Fisher matrix, i.e. C = F and the root mean square error of the 
parameters can be evaluated from the diagonal elements as 


\/((AA)2) = v^. 


(17) 


When analysing the synergy effect of combining the measurements from two detectors, we take the sum of the Fisher 
matrices corresponding to them and then take the inverse of the matrix thus obtained: 

Ecomb =ri+r2, (18) 

Ccomb = Tcomb- (1^) 
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Studying the change in estimates by varying design sensitivity of DECIGO 

We also study the synergy effect on the error estimate of various parameters by varying the sensitivities of the 
future space mission DECIGO (and its precursor pre-DECIGO), since its design sensitivity is still to be determined 
depending on the scientific gain. Here, for brevity, we simply call it DECIGO. In both better and worse directions, 
the sensitivities are changed by scaling the fiducial noise curve uniformly over all frequencies. Namely, the scaled 
DECIGO noise curves are obtained as 

S^(^f)scaled^l^SM)DECIGO^ ( 20 ) 

with a constant 1C. 


III. RESULTS 

In this section we explore the expected errors in the parameter estimation when we combine outputs of ground and 
space based detectors. We discuss results on the various parameters in the PN and the ppE waveforms for coalescing 
binary systems. As mentioned earlier we consider two NS-BH binary systems: a 1.4 Mgoiar + 10 Mgoiar binary and a 
1.4 Msolar + 100 Msolar binary with the distance fixed at 200 Mpc. We fix the lower frequency limits at 0.1 Hz for 
the space based GW detectors. The space based detectors are expected to be sensitive at frequencies much below this 
cut off (eLISA is expected to be sensitive at as low as 0.1 mHz), but we chose this cut off since the evolution of the 
binary coalescence is really slow at lower frequencies. For the mass configurations we have considered in this work, 
it will take many years for the binaries to span the lower (10“® - 0.1) frequency range and we would not be able to 
observe this continuously. The evolution is much faster in the mid-high frequency range and there is a chance we can 
observe the same events in both space and ground based detectors. For advanced-LIGO we set the lower frequency 
limit at 10 Hz. 


A. Parameter estimation - PN expansion 


Here we present the error estimates on the chirp mass, the symmetric mass ratio, and the phase and time at 
coalescence using the method described in the previous section. In table I and H we tabulate these errors for both 
independent and combined measurements. We note that there is some gain in combining space-ground measurements 
for all parameters. In case of the chirp mass, maximum gain is obtained for the detector pair Advanced-LIGO - eLISA, 
while for all the other pairs one of the GW detector completely dominates the error estimate. Similarly for the time 
of coalescence maximum gain is obtained for the detector pair ET-DECIGO, and for the other pairs the estimates 
are dominated by one of the two detectors. For the symmetric mass ratio and the phase of coalescence there is no 
substantial gain in combining the space-ground measurements since one of the detector always dominates the results. 
We further analyse the variation in the error estimates on the parameters with the variation of the DEGIGO sensitivity 
according to equation (20). As a reference, we introduce a quantity A0' for a parameter 0, which corresponds to the 
error estimate from the combined measurement if that single parameter was being estimated (equation (21)). In such 
a case, the joint probability for the combined estimate of some parameter say 9 is given by 

P{e) (X X 

^ g-(Ae)V(2(Ae')")^ 


where A0i and A 02 are the error estimates from the first and the second measurements, respectively, and 


M' = 


A6>iA6i2 

^A02 + • 


( 21 ) 


These reference values are, off course, different from the correct estimate A0 (combined) = based on Ccomb 

given in equation (19). The correct estimates of errors normalized by these reference values are shown in figure (2) 
and (3). From our study on the synergy between (scaled-)DECIGO and ET, we can infer the following about the PN 
template parameters: for the chirp mass, there is a wide range of sensitivities in the pre-DECIGO phase (sensitivities 
that are 10 - 10^ order of magnitude worse than DEGIGO) where there is substantial gain in combining the space- 
ground measurements, for the symmetric mass ratio and the phase at coalescence maximum gain is obtained when 
the pre-DECIGO sensitivity is almost 10^ order of magnitude worse than DECIGO, and for the time at coalescence 
maximum gain is achieved in the post-DEGIGO phase where the sensitivity is almost an order of magnitude better 


than DECIGO. 
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with scaled-DECIGO noise curve:1.4 Msoia^+lOO M^aiar 



(a) 


with scaled-DECIGO noise curve:1.4 Msoiar+100 Msoiar 



^^DECIGo/ 

(b) 


FIG. 2: Variation in the error estimates on the chirp mass and the symmetric mass ratio (a), and the time and 
phase at coalescence (b), obtained by varying DECIGO sensitivities. A0comb = M-jvjiclvc^ is obtained as 
explained in section IID, and A0' is as given in equation (21). The ratio A0comb/A0' is plotted against A0 as 
obtained from the scaled DECIGO measurements. The dotted curves correspond to better sensitivity than DECIGO 


with scaled DECIGO noise curves (1.4 Mg+100 M^j 



.^(f)Scaled/^^^(f)DECIGO 

(a) 


with scaled DECIGO noise curves (1.4 Mg+lOO Mg) 



(b) 


EIG. 3: Variation in the error estimates on the chirp mass and the symmetric mass ratio (a), and the time and 
phase at coalescence (b), with varying DECIGO sensitivities. A0comb = -Mfvjtclpc) is obtained as explained in 
section IID, and A9' is as given in equation (21). Here too, the dotted curves correspond to better sensitivity than 

DECIGO. 


B. ppE expansion 

We present the results of the analysis for IPN, 2PN and 3PN modifications to GR in this section. Here IPN, 
2PN and 3PN modihcations mean that we hx b = —1, b = —1/3 and 6 = 1/3 in equation (8), respectively. The 
IPN correction corresponds to the traditional Massive graviton theory ([11]), and the 2PN correction corresponds 
to Quadratic curvature theory ([11]). With distance fixed at 200 Mpc the results are shown in table HI and IV for 
IPN, table V and VI for 2PN and table VH and VHI for 3PN. In all cases we find some gain to the constraint on 
these modihcations to GR when combining the measurements of the space and ground based detectors. In order to 
clarify the synergy effect, we present the plots similar to 2 and 3 for the parameters of the magnitude of modihcation 










with scaled DECIGO noise curves (1.4 Mg+lOO Mg) 


with scaled DECIGO noise curves (1.4 Mg+100 Mg) 




(a) (b) 

FIG. 4: Variation in the error estimates on the GR modification parameter at different PN orders. A/3 (combined) is 
obtained as explained in section IID, and A/3' is as given in equation (21). The dashed part corresponds to the 

noise curves with better sensitivity than DEGIGO 


with scaled DECIGO noise curves (1.4 Mg+100 Mg) 


with scaled DECIGO noise curves (1.4 Mg+100 Mg) 




FIG. 5: Variation in the error estimates on the GR modification parameter with varying pre-DEGIGO sensitivity. 
A/3 (combined) is obtained as explained in section IID, and A/3' is as given in equation (21). The dotted curves 

correspond to better sensitivity than DECIGO. 


from GR. As before, in figure (4) we plot the error estimate based on Ccomb given in equation (19), IS.j3, relative 
to the reference value A/3' defined by equation (21), for combined measurements from scaled-DECIGO and ET. For 
very low DECIGO sensitivity the error estimates are dominated by ET, while for very high DECIGO sensitivity 
they are dominated by DECIGO. In both limiting cases we do not have any merit of combining the space-ground 
measurements. However, in the middle range we have a sweet spot where the gain becomes significantly large. In 
this plot (figure (4)) the horizontal axis is A,5decigo/A/3et- For all PN order modifications, the sweet spot appears 
around A/3decigo/^/3et = 1; where both detectors equally contribute to constrain the parameter /3. In figure (4b) 
we plot I^Pcombined- In the region with A/3decigo/A/3et "C 1, the error estimates are dominated by the DECIGO 
sensitivity, and hence they show linear dependence. In the region with A/3decigo / A,0et ^ 1, the error estimates 
are dominated by the ET sensitivity, and hence the curves become almost flat. At around A/3decigo/A/3et = 1, the 
slope of the curves changes. This change is not monotonic, and corresponds to the sweet spot in the figure (4a) where 
the true joint estimates are better than the reference estimates given by A/3'. 













9 


We also show these error estimates as functions of the sensitivity level of the scaled-DECIGO in figure (5). The 
higher order post-Newtonian terms tend to be better determined by the ground based detectors which have better 
sensitivities at higher frequencies. Hence, for higher PN order terms, the synergy is maximum (i.e. the space and 
ground detectors contribute almost equally) for higher sensitivity of scaled-DECIGO. The positions of the sweet spots 
in figure (3) can be understood in the same way by noticing that the leading order frequency dependences of the chirp 
mass, the symmetric mass ratio, and the phase and time at coalescence are OPN, IPN, 2.5PN and 4PN, respectively. 


IV. CONCLUSIONS 

In this paper we have assessed the expected synergy effects between ground and space based detectors on the 
parameter estimation of coalescing binary systems. Our aim is to demonstrate that the advantage of having a GW 
antenna which is sensitive at low frequencies, is larger than we naively expect. Larger gain in the error estimate of 
the GW waveform parameters is obtained when the constraints from individual ground or space detectors are almost 
equal. In case of the ppE parameters that characterize deviations from GR, some gain is always obtained irrespective 
of the level of sensitivity of the GW antenna. For lower PN order modifications some gain can be obtained even 
if the sensitivity of the space detector is degraded. The same argument also applies to the estimated errors in the 
extraction of binary parameters contained in the standard PN templates. In the present paper, neglecting the spins, 
we considered the chirp mass, the symmetric mass ratio, and the phase and time at coalescence, whose leading order 
frequency dependences on the GW phase are OPN, IPN, 2.5PN and 4PN, respectively. We found a sweet spot in the 
combined estimates of the chirp mass and the symmetric mass ratio in the pre-DEGIGO phase (worse sensitivity than 
DEGIGO). In this region, the combined error estimates are better than individual detector estimates. Substantial gain 
was also obtained in the joint measurement of the time and phase at coalescence. Here too the sweet spot was in the 
pre-DEGIGO phase for the phase at coalescence, but for the time at coalescence, maximum gain was obtained with 
post-DEGIGO sensitivity. As expected, for higher PN order parameters, synergy is obtained at higher sensitivities of 
the space-borne GW detector. 

The sensitivity of space-borne GW antenna such as pre-DEGIGO, which is the precursor mission for DEGIGO, is 
still to be determined by considering various aspects, including the scientific gain. The study of the synergy effects 
reported in the present paper provides valuable information and can be taken into account to make such decisions. 
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Appendix: Tabels 


TABLE I: Errors estimates for PN waveform parameters for individual detector measurements. These are calculated 
for NS-BH binaries with masses 1.4 Mq-\-10 Mq and 1.4 Mq-1-100 Mq. The distance is fixed to 200 Mpc. 


Binary mass 

Detector 

AC 

A</>c 

AM/M(%) 

Av/v{%) 

1.4M0-hlOM0 

DECIGO 

9.16 X 10“® 

9.70 X 10-^ 

9.57 X 10-® 

2.60 X 10"^ 

1.4M0-hlOM0 

eLISA 

1.30 X 10^ 

3.87 X 10^ 

8.16 X 10"® 

2.80 X 10^ 

1.4M0-hlOM0 

advLIGO 

3.22 X 10“'‘ 

3.26 X 10"^ 

1.62 X 10"^ 

4.52 X 10"^ 

1.4M0-hlOM0 

ET 

1.51 X 10“® 

1.64 X 10"^ 

1.35 X 10"'‘ 

1.31 X 10"^ 

1.4M0-hlOOM0 

DECIGO 

6.52 X 10“® 

2.43 X 10"® 

7.36 X 10"® 

9.36 X 10"® 

1.4M0-hlOOM0 

eLISA 

1.08 X 10^ 

2.07 X 10^ 

6.20 X 10"^ 

6.73 X 10^ 

1.4M0-hlOOM0 

advLIGO 

2.09 X 10-^ 

1.49 X 10"^ 

3.43 X 10"^ 

3.53 X 10"^ 

1.4M0-hlOOM0 

ET 

1.02 X 10“® 

6.02 X 10"® 

3.00 X 10"^ 

1.23 X 10"® 


TABLE II: Errors estimates for PN waveform parameters for combined detector measurements. These are calculated 
for NS-BH binaries with masses 1.4 Mq-\-10 Mq and 1.4 Mq-1-100 Mq. The distance is fixed to 200 Mpc. 


Binary mass 

Detector 

Ate 

<1 

AMIM{%) 

Aulu{%) 

1.4M0-I-1OM0 

ET-hDECIGO 

6.04 X 10"® 

8.98 X 10"^ 

9.15 X 10“® 

2.46 X lO""* 

1.4M0-I-1OM0 

ET-heLISA 

1.47 X 10"® 

1.58 X 10"^ 

1.03 X 10“'‘ 

1.21 X 10"® 

1.4M0-I-1OM0 

advLIGO-l-DEGIGO 

4.95 X 10"® 

9.22 X 10"^ 

9.28 X 10“® 

2.50 X lO""* 

1.4M0-I-1OM0 

advLIGO-heLISA 

2.36 X 10"'‘ 

2.55 X 10"^ 

1.79 X 10-^ 

2.41 X 10"® 

1.4M0-hlOOM0 

ET-hDECIGO 

3.31 X 10"® 

1.20 X 10"® 

3.92 X 10“® 

5.24 X 10"® 

1.4M0-hlOOM0 

ET-heLISA 

9.84 X 10"® 

5.73 X 10"® 

2.28 X 10-^ 

1.14 X 10"® 

1.4M0-hlOOM0 

advLIGO+DEGIGO 

2.56 X 10"® 

1.56 X 10"® 

4.90 X 10“® 

6.52 X 10"® 

1.4M0-hlOOM0 

advLIGO-heLISA 

1.74 X lO""' 

1.36 X 10"® 

3.53 X 10“^ 

2.22 X 10"® 





TABLE III: Errors estimates for ppE waveform parameters for individual detector measurements, when the 
modification to GR appears at the first PN order. These are calculated for NS-BH binaries with masses 

1.4 Mq+100 Mq. The distance is fixed to 200 Mpc. 


11 


Binary mass 

Detector 

Ate 

Abe 

AA4/A4(%) 

Avlv{%) 

A/3 

1.4M0-MOM0 

DECIGO 

1.56 X 

9.37x10"® 1.85x10"^ 

9.78 X 10-^ 

2.48 X 10"'^ 

1.4M0-MOM0 

eLISA 

2.59 X 10^ 

4.91 X 10® 

2.98 X 10“® 

3.00 X 10® 

1.15 X 10"® 

1.4M0-MOM0 

aLIGO 

7.16 X 10"^ 

7.09 X 10"® 9.13 X 10“® 

2.33 

8.48 X 10"® 

1.4M0-MOM0 

ET 

2.33 X 10"® 

1.75 X 10"® 4.64 X 10“'‘ 

4.44 X 10"® 

8.86 X 10"® 

1.4M0-hlOOM0 

DECIGO 

7.56 X 10"® 

3.55 X 10"® 3.76 x 10"^ 

2.03 X lO""* 

1.69 X 10"'^ 

1.4M0-hlOOM0 

eLISA 

1.22 X 10^ 

2.12 X 10® 

9.18 X 10“® 

6.74 X 10® 

1.26 X 10"® 

1.4M0-hlOOM0 

aLIGO 

7.18 X 10"'‘ 

1.93 

1.32 X 10“® 

1.16 X 10"® 

5.46 X 10"® 

1.4M0-hlOOM0 

ET 

2.01 X 10"® 

3.42 X 10"® 7.98 x 10-^ 

3.36 X 10"® 

7.89 X 10"® 


TABLE IV: Errors estimates for ppE waveform parameters for combined detector measurements, when the 
modification to GR appears at the first PN order. These are calculated for NS-BH binaries with masses 

1.4 Mq + 100 Mq. The distance is fixed to 200 Mpc. 


Binary mass 

Detector 

Ate 

Abe 

AA4/A4(%) 

Aulu{%) 

A/3 

1.4M0-hlOM0 

ET-hDECIGO 

7.43 X 10"® 

4.54 X 10"® 

1.33 X 10"'^ 

5.84 X 10"® 

1.26 X lO"®" 

1.4M0-hlOM0 

ET-heLISA 

1.91 X 10“® 

1.72 X 10"® 

1.82 X 10"® 

2.86 X 10"® 

4.50 X 10“® 

1.4M0-hlOM0 

aLIGO+DECIGO 

5.71 X 10"® 

6.12 X 10"® 

1.55 X 10"'^ 

7.32 X 10"® 

1.67 X lO"®" 

1.4M0-hlOM0 

aLIGO-beLISA 

3.3 X lO""® 

2.60 X 10"® 

1.10 X 10"® 

6.10 X 10"® 

1.13 X 10“® 

1.4M0-MOOM0 

ET-hDECIGO 

3.87 X 10"® 

1.47 X 10"® 

2.53 X 10"'^ 

1.09 X 10"® 

1.30 X lO"®" 

1.4M0-MOOM0 

ET-heLISA 

1.59 X 10“® 

2.23 X 10"® 

3.66 X 10"® 

2.50 X 10“® 

4.72 X 10“® 

1.4M0-MOOM0 

aLIGO+DECIGO 

2.67 X 10"® 

2.15 X 10"® 

3.09 X 10"'^ 

1.46 X 10"® 

1.51 X lO"®" 

1.4M0-MOOM0 

aLIGO+eLISA 

3.01 X 10“® 

4.95 X 10"® 

1.85 X 10"® 

5.14 X 10"® 

1.21 X 10“® 


TABLE V: Errors estimates for ppE waveform parameters for individual detector measurements, when the 
modihcation to GR appears at the second PN order. These are calculated for NS-BH binaries with masses 

1.4 Mq-1-100 Mq. The distance is fixed to 200 Mpc. 


Binary mass 

Detector 

Ate 

Abe 

AM/M{%) 

Aulu(%) 

A/3 

1 . 4 M 0 + 1 OM 0 

DECIGO 

1.63 X 10"® 

1.27 X 10"® 

1.66 X lO"®" 

5.62 X lO""® 3.51 X 10"® 

1 . 4 M 0 + 1 OM 0 

eLISA 

2.76 X 10® 

7.27 X 10® 

2.47 X 10“® 

1.17 X 10® 

1.84 X 10® 

1 . 4 M 0 + 1 OM 0 

aLIGO 

3.95 X 10"® 

8.21 X 10® 

7.15 X 10“® 

4.21 X 10® 

1.60 X 10® 

1 . 4 M 0 + 1 OM 0 

ET 

1.00 X 10"® 

1.90 

1.21 X 10“® 

9.95 X 10"® 3.72 x 10"® 

1 . 4 M 0 + 1 OOM 0 

DECIGO 

7.45 X 10"® 

2.53 X 10“® 

4.40 X lO"®" 

3.57 X lO""® 3.88 x 10~® 

1 . 4 M 0 + 1 OOM 0 

eLISA 

1.22 X 10® 

2.20 X 10® 

6.26 X 10“® 

6.76 X 10® 2.75 x 10® 

1 . 4 M 0 + 1 OOM 0 

aLIGO 

1.18 X 10"® 

7.02 

6.25 X 10“® 

3.19 X 10"® 9.05 X 10"® 

1 . 4 M 0 + 1 OOM 0 

ET 

4.42 X 10"® 

2.22 X 10"® 

4.08 X lO"'® 

1.14 X 10"® 2.78 X 10"® 





TABLE VI: Errors estimates for ppE waveform parameters for combined detector measurements, when the 
modification to GR appears at the second PN order. These are calculated for NS-BH binaries with masses 

1.4 Mq+100 Mq. The distance is fixed to 200 Mpc. 
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Binary mass 

Detector 

Ate 

Abe 

AA4/A4(%) 

Aulu{%) 

A/3 

1.4M0-hlOM0 

ET-hDECIGO 

7.64 X 10“® 

5.80 X 10"® 

1.21 X 10"'^ 

3.76 X 10"^ 

1.66 X 10“'‘ 

1.4M0-hlOM0 

ET-heLISA 

3.27 X 10“® 

4.88 X 10"^ 

2.50 X 10"^ 

2.62 X 10"^ 

9.65 X 10-2 

1.4M0-hlOM0 

aLIGO+DECIGO 

5.77 X 10“® 

8.09 X 10"^ 

1.41 X 10"'^ 

4.54 X 10"^ 

2.28 X 10-^ 

1.4M0-hlOM0 

aLIGO-^eLISA 

5.20 X 10-^ 

8.23 

1.83 X 10"® 

4.37 

1.61 

1.4M0-MOOM0 

ET-hDECIGO 

4.15 X 10“® 

2.04 X 10"^ 

3.13 X 10"'^ 

2.36 X 10"^ 

3.06 X 10-^ 

1.4M0-MOOM0 

ET-heLISA 

4.03 X 10“® 

1.96 X 10"^ 

2.73 X 10"'‘ 

1.03 X 10"^ 

2.44 X 10“2 

1.4M0-MOOM0 

aLIGO+DEGIGO 

2.67 X 10“® 

2.12 X 10"® 

3.74 X 10"'^ 

2.92 X 10"^ 

3.52 X 10“'‘ 

1.4M0-MOOM0 

aLIGO-^eLISA 

7.57 X lO-'^ 

3.91 

5.42 X 10"'‘ 

1.98 X 10"^ 

4.97 X 10“^ 


TABLE VII: Errors estimates for ppE waveform parameters for individual detector measurements, when the 
modification to GR appears at the third PN order. These are calculated for NS-BH binaries with masses 

1.4 Mq+100 Mq. The distance is fixed to 200 Mpc. 


Binary mass 

Detector 

Ate 

<1 

AA4/A4(%) 

Aulu{%) 

A/3 

1 . 4 M 0 + 1 OM 0 

DEGIGO 

2.29 X lO""' 

6.43 X 10"® 1.29 X 10"^ 

3.86 X 10-^ 2.48 x 10-^ 

1 . 4 M 0 -MOM 0 

eLISA 

3.98 X 102 

4.68 X lO^ 

1.65 X 10-2 

6.59 X 10^ 

2.08 X 10® 

1 . 4 M 0 -MOM 0 

aLIGO 

1.19 X 10-® 

7.69 

2.93 X 10-2 

2.44 

4.67 X 10^ 

1 . 4 M 0 + 1 OM 0 

ET 

4.80 X 10"® 

2.15 X 10"^ 1.52 X 10“'‘ 

5.08 X 10-2 

1.38 

1 . 4 M 0 -MOOM 0 

DEGIGO 

1.8 X 10"^ 

3.24 X 10-2 6.83 x 10-^ 

6.75 X IQ-"* 4.02 x 10-^ 

1 . 4 M 0 -MOOM 0 

eLISA 

1.38 X 102 

2.50 X 10^ 

6.51 X 10-2 

7.02 X 10^ 

3.07 X 10^ 

1 . 4 M 0 -MOOM 0 

aLIGO 

2.09 X lO""' 

2.96 

7.44 X 10-2 

3.79 X 10-^ 

2.97 X 10^ 

1 . 4 M 0 -MOOM 0 

ET 

1.28 X 10"® 

7.81 X 10-2 4.12 X 10-^ 

7.53 X 10-® 7.17 X 10-^ 


TABLE VIII: Errors estimates for ppE waveform parameters for combined detector measurements, when the 
modification to GR appears at the third PN order. These are calculated for NS-BH binaries with masses 

1.4 Mq- flOO Mq. The distance is fixed to 200 Mpc. 


Binary mass 

Detector 

Ate 

Abe 

AA4/A4(%) 

Aulu{%) 

A/3 

1 . 4 M 0 -MOM 0 

ET-bDECIGO 

1.15 X 10-® 

2.15 X 10-® 

9.73 X 10-® 

2.69 X 10-^ 

6.62 X 10-2 

1 . 4 M 0 -MOM 0 

ET-beLISA 

4.60 X 10-® 

1.98 X 10-^ 

1.07 X 10-'‘ 

4.53 X 10-2 

1.28 

1 . 4 M 0 -MOM 0 

aLIGO+DEGIGO 

6.72 X 10-® 

3.70 X 10-® 

1.13 X 10-'^ 

3.27 X 10-^ 

1.3 X 10-^ 

1 . 4 M 0 -MOM 0 

aLIGO+eLISA 

8.41 X 10-^ 

3.99 

3.85 X 10-^ 

1.00 

2.52 X 10^ 

1 . 4 M 0 -MOOM 0 

ET+DECIGO 

7.10 X 10-® 

5.65 X 10-® 

1.30 X 10-'^ 

1.44 X 10-^ 

5.53 X 10-2 

1 . 4 M 0 -MOOM 0 

ET+eLISA 

1.27 X 10-® 

6.75 X 10-2 

2.67 X lO-"' 

6.2 X 10-® 

6.11 X 10-^ 

1.4Mq-MOOM0 

aLIGO+DEGIGO 

4.31 X 10-® 

1.77 X 10-2 

3.87 X 10-'^ 

3.94 X 10-^ 

2.00 X 10-^ 

1 . 4 M 0 -MOOM 0 

aLIGO+eLISA 

2.03 X lO-'^ 

1.41 

3.58 X 10-^ 

1.50 X 10-^ 

1.36 X 10^ 





